Application of Information Theory in Nuclear Liquid Gas Phase Transition 
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Information entropy and Zipf 's law in the field of information theory have been used for studying 
the disassembly of nuclei in the framework of the isospin dependent lattice gas model and molecular 
dynamical model. We found that the information entropy in the event space is maximum at the 
phase transition point and the mass of the cluster show exactly inversely to its rank, i.e. Zipf's 
law appears. Both novel criteria are useful in searching the nuclear liquid gas phase transition 
experimentally and theoretically. 
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Hot nuclei can be formed in energetic heavy ion colli- 
sions (HIC) and may be highly excited. They deexcite 
by different decay modes, such as multifragmentation. 
Experimentally, this kind of multifragment emission was 
observed to evolve with beam energy (excitation energy, 
or nuclear temperature, ...). Multiplicity, Ni^f, of inter- 
mediate mass fragment (IMF) rises with the beam energy, 
reaches a maximum, and finally falls to lower value. This 
phenomenon of the rise and fall of Nimf may be related 
to the liquid gas phase transition in nuclear matter [Q. 
The onset of multifragmentation probably indicates the 
coexistence of liquid and gas phases. The mass distribu- 
tion of IMF distribution can be expressed as power law 
with parameter r. The minimum of t, Tmim occurs when 
the liquid gas phase transition takes place However, 
Tmin Can also reveal at supercritical densities along the 
Kertesz line B and at some subcritical densities at lower 
temperature IJ] . So it is not possible to determine the 
phase transition only from A^im/ find Tmin- 

On the other hand, experimentalists measured the nu- 
clear caloric curves, i.e. the relationship between nuclear 
temperature and the excitation energy. He-Li isotopic 
temperature from Albergo thermometer for projectile- 
like Au spectators seems to exhibit a temperature plateau 
in the excitation energy range of 3 to 10 MeV/u g. This 
plateau was taken as an indication for a first order nu- 
clear liquid gas phase transition. However, due to the 
changing mass of Au spectators with excitation energy 
and the side-feeding effect to measured He-Li isotopic 
temperature, this conclusion is questionable |Q. Nu- 
clear caloric curves were also surveyed by several groups 
1^. However, unfortunately, the sharp signature of liq- 
uid gas phase transition in macroscopic systems may be 
smoothed and blurred due to the small numbers of nu- 
cleons in nuclei, and/or the difficulty to perform a di- 
rect comparison between the measured " apparent" tem- 
perature and the "real" temperature interferes obtaining 
the real nuclear caloric curve. These factors hamper the 
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reaching of a definite conclusion on liquid gas phase tran- 
sition in nuclei. 

The extraction of critical exponents and the study of 
critical behavior in finite-size systems were atteniptcd in 
and were followed by controversial debates |jl^. In 
this context, it is necessary and meaningful to search for 
some novel signatures to characterize the nuclear liquid 
gas phase transition in order to guide the experimental 
analysis and theoretical predictions. 

In this Letter, we will introduce information entropy 
]ll[ , H, and Zipf's law |Q into the diagnosis of nuclear 
liquid gas phase transition. The information entropy was 
defined by Shannon in information theory. Originally, it 
measures the " amount of information" which is contained 
in messages sent along a transmission line. It can be 
expressed as follows, 

H=-Y,PMP^)■ (1) 

i 

where pi is a normalized probability, and X^iP* ~ 1- 
Jaynes proposed that a very general technique for dis- 
covering the least biased distribution of the pi consists 
in the maximization of the Shannon H entropy, subject 
to whatever constraints on pi are appropriate to the par- 
ticular situation. The maximization of H was thus put 
forward as a general principle of statistical inference - 
one which could be applied to a wide variety of prob- 
lems in economics, engineering and many other fields, 
such as quantum phenomena pi| . In high energy hadron 
collisions, multiparticle production proceeds on the max- 
imum stochasticity, i.e. they should obey the maximum 
entropy principle. This kind of stochasticity can be also 
quantified via the information entropy which has been 
shown to be a good tool to measure chaoticity in hadron 
decaying branching process ||l^. In different physical 
conditions, information entropy can be expressed with 
different stochastic variables. In this work on HIC, we 
define pi as the event probability of having " z" particles 
produced, i.e. {pi} is the normalized probability distri- 
bution of total multiplicity, the sum is taken over whole 
{pi} distribution. This emphasis is on the event space 
rather than the phase space. As shown below this kind 
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of information entropy |TJ] can be taken as a method to 
determine nuclear liquid gas phase transition. 

Zipf's law [|l2j has been known as a statistical phe- 
nomenon concerning the relation between English words 
and their frequency in literature in the field of linguistics. 
The law states that, when we list the words in the order 
of decreasing population, the frequency of a word is in- 
versely proportional to its rank |12| . This relation was 
found not only in linguistics but also in other fields of 
sciences. For instance, the law appeared in distributions 
of populations in cities, distributions of income of cor- 
porations, distributions of areas of lakes and cluster-size 
distribution in percolation process p5[ |. In this Letter, 
Zipf's law will be tested for the fragment mass distribu- 
tion and it is evidenced to be a factor to characterize the 
phase transition. 

The tools we will use here are the isospin dependent 
lattice gas model (LGM) and molecular dynamical model 
(MD). The lattice gas model was developed to describe 
the liquid-gas phase transition for atomic system by Lee 
and Yang The same model has already been ap- 

plied to nuclear physics for isospin symmetrical systems 
in the grandcanonical ensembl e w ith a sampling of 
the canonical ensemble U El fsTll |20l 1211, 



also for isospin asymmetrical nuclear matter in the mean 
field approximation [^ . In addition, a classical molecu- 
lar dynamical model is used to compare its results with 
the results of lattice gas model. Here we will make a brief 
description for the models. 

In the lattice gas model, A N + Z) nucleons with 
an occupation number s which is defined s = 1 (-1) for 
a proton (neutron) or s = for a vacancy, are placed on 



the L sites of lattice. Nucleons in the nearest neighbor- 
ing sites interact with an energy es-s^- The hamiltonian 

is written as E J2f=i ^-J2i<j (^siSjS^Sj. The interac- 
tion constant es^sj is chosen to be isospin dependent and 
be fixed to reproduce the binding energy of the nuclei 



= O.MeF, 



= -b.SSMeV. 



(2) 



A three-dimension cubic lattice with L sites is used. The 
freeze-out density of disassembling system is assumed to 
be pf — j;Po, where po is the normal nuclear density. 
The disassembly of the system is to be calculated at p/, 
beyond which nucleons are too far apart to interact. Nu- 
cleons are put into lattice by Monte Carlo Metropolis 
sampling. Once the nucleons have been placed we also 
ascribe to each of them a momentum by Monte Carlo 
samplings of Maxwell-Boltzmann distribution. 

Once this is done the LGM immediately gives the clus- 
ter distribution using the rule that two nucleons are part 



and of the same cluster if P^/2/x — es;s -SiSj < 0. This method 



is similar to the Coniglio-Klein's prescription |^ in con- 
densed matter physics and was shown to be valid in LGM 
TTI. To calculate clusters using MD we propa- 
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gate the particles from the initial configuration for a long 
time under the influence of the chosen force. The form of 
the force is chosen to be also isospin dependent in order 
to compare with the results of LGM. The potential for 
unlike-nucleons is 



Unp{r){r/ro <a) = C [Biro/rf - {ro/ry]exp{ 



1 



Vnp{r){r/ro > a) 



0. 



(r/ro) - a 



(3) 



r 



where ro = 1.842 fm is the distance between the centers 
of two adjacent cubes. The parameters of the potentials 
are p ^ 2, q = 1, a = 1.3, B = 0.924, and C = 1966 
MeV. With these parameters the potential is minimum 
at ro with the value -5.33 MeV, is zero when the nucleons 
are more than 1.3ro apart and becomes strongly repulsive 
when r is significantly less than tq. The potential for like- 
nucleons is written as 



'pp(r)[r < ro) 



Vppir){-. 
Vppir){r > 



= Vnpir) ~ Vnpiro), 



ro) = 0. 



(4) 



The system evolves with the above potential. At asymp- 
totic times the clusters are easily recognized. Observables 
based on the cluster distribution in the both models can 
now be compared. In the case of proton-proton inter- 
actions, the Coulomb interaction can also be added sep- 



arately and it can be compared with the case without 
Coulomb effects. 

In this Letter we choose the medium size nuclei -'^^^Xe 
as an example. In most cases, pf is chosen to be 0.38 
Po, since the experimental data can be best fitted by p/ 
between 0.3po and OApo in previous LGM calculations 
Ip^ , |25[| , which corresponds to 7^ cubic lattice. In addi- 
tion, O.lSpo, corresponding to 9"^ cubic lattice and 0.60po, 
corresponding to 6^^ cubic lattice of p/ are also taken to 
compare and check the results with different pf values 
in the LGM case. 1000 events were simulated for each 
combination of T and p f which ensures enough statistics 
for results. 

In order to check the phase transition behavior in LGM 
and MD, we will first show the results of some physical 
observables, namely the effective power-law parameter, r. 
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FIG. 1: The effective power-law parameter, r, the second moment 
of the cluster distribution, S2, and the multiplicity of intermediate 
mass fragments, Ni^f as a function of temperature for the disas- 
sembly of ^^^Xe. Left panel is the LGM results with different pf 
and right panel is the comparison of MD to LGM with 0.38po- The 
symbols are illustrated in figure. 

the second moment of the cluster distribution, 5*2, and 
the multiphcity of intermediate mass fragments, Nimf, 
for the disassembly of ^^^Xe: fig.l. These observables 
were shown to be good indicators of a liquid gas phase 
transition, as shown in Ref. pO| . The extreme values 
of r, Nimf and S2 occur at the same temperature, in- 
dicating of the onset of the phase transition, for each 
calculation case. For the LGM case, the phase transi- 
tion temperature increases with the freeze-out density; 
for the MD case, a slight small transition temperature 
is obtained when Coulomb force is ignored. It becomes 
much lower in case of Coulomb interaction due its long 
range repulsion. Similar phenomenon has been explored 
in a previous study | pO| . However, the aim of this Letter 
is to take the above transition temperatures as references 
to search for novel signatures of liquid gas phase transi- 
tion. 



Fig. 2 shows the information entropy for disassembly of 
Xe. The information entropy exhibits a rise and fall with 
temperature, which is similar to the behaviors of Nimf 



FIG. 2: The same as Fig.l, but for the information entropy H. 

and 5*2. The temperatures extracted from the peak val- 
ues of H are consistent with the transition temperatures 
in Fig.l, indicating that information entropy ought to 
be a good diagnosis of phase transition. Physically, the 
maximum of H reflects the largest fluctuation of the mul- 
tiplicity probability distribution in the phase transition 
point. In this case it is the most difficult to predict how 
many clusters will be produced in each event, i.e. the 
disorder (entropy) of information is the largest. Gen- 
erally speaking, the larger the dispersal of multiplicity 
probability distribution, the higher the information en- 
tropy and then the disorder of system in the event topol- 
ogy. One should make a careful distinction between this 
information entropy, on the one hand, and the original 
thermodynamic entropy, on the other hand |l^ . The 
latter generally illustrates the heat disorder in momen- 
tum space rather than event space and it always increases 
with temperature. 

Next we will present the results for testing Zipf's law 
in the fragment distribution. The law states that the re- 
lation between the sizes and their ranks is described by 
An — c/n (n=l, 2, 3, ...), where c is a constant and An 
is the mass of rank n in a mass list when we arrange 
the clusters in the order of decreasing size. In the cal- 
culations, we averaged the masses for each rank in mass 
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lists of the events. Then we examined the relation be- 
tween the masses An and their ranks n with the fit of 
An oc n^^ in the range of 1 < n < 10, where A is the 
slope parameter. The upper panel of Fig. 3 summarizes 
such parameter A as a function of temperature for LGM 
and MD cases. Clearly, the value of A decreases with tem- 
perature, indicating that the difference of mass between 
the different fragment ranks is becoming smaller. When 
A 1, the Zipf's law is satisfied: An oc n^^. The tem- 
peratures having Zipf's law for all calculations in Fig. 3 
are also consistent with the respective transition temper- 
atures extracted from the extreme values of some observ- 
ables in Fig.l and 2. Therefore Zipf's law is also a good 
signal to phase transition. From the statistical point of 
view, Zipf's law is related to the critical behavior or self- 
organized criticality ^ , which may be a special state 
with the maximum information or least effort. 

In order to further illustrate that Zipf's law exists at 
the phase transition point most probably, we directly fit 
the rank-classified cluster distribution with Zipf's law 
and extract the truth of the hypothesis: test. The 
lower panel of Fig. 3 demonstrates the /ndf for the An 
- n relations at different T for different calculation cases. 
As expected, there are minima of /ndf around the re- 
spective transition temperature, which further support 
Zipf's law of the fragment distribution indicates a phase 
transition. All calculations give the same conclusions as 
above. 



its rank, i.e. Zipf's law appears. Even thought both 
criteria are still phenomenological, we believe that they 
are simple and practicable tools to diagnose the nuclear 
liquid gas phase transition in experiments and theories. 
We are waiting for some data analysis on the information 
entropy and Zipf's law in near future. 
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In addition, we also investigated larger systems, such 
as A = 274, 500 and 830, in the LGM case to see if the 
system behaves as expected in ^^^Xe. The results show 
that the maximum of information entropy and Zipf's law 
behavior still remain at the same phase transition tem- 
perature as the one extracted from the extreme values of 
T, Nirnf and 52- It illustrates that both novel criteria are 
suitable as signals of phase transition in larger A limit. 

In conclusion, the information entropy and Zipf's law 
in the field of information theory are introduced, for the 
first time, into the study of the liquid gas phase transi- 
tion of nuclei in the framework of the isospin dependent 
lattice gas model and molecular dynamical model in p-T 
plane. At the point of phase transition, the information 
entropy of multiplicity distribution is maximum, which 
indicates that the system at this time has the largest 
fiuctuation/stochasticity/chaoticity in the event space. 
Meanwhile, the cluster mass show exactly inversely to 



FIG. 3: The slope parameter A of the relation An to n (upper 
panel) and the x^/nd/ with the fit of Zipf's law (lower panel) as 
a function of temperature. Left panel is the LGM results with 
different pf and right panel is the comparison of MD to LGM with 
0.38po- The symbols are the same as Fig.l. 
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